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Abstract 

Based on the quantum hydrodynamic model, a new relationship between the 
electrostatic-potential and the electron-density in the ultra-dense plasma is derived. 
Propagation of arbitrary amplitude nonlinear ion waves is, then, investigated in a com- 
pletely degenerate dense dusty electron-ion plasma, using this new energy relation, 
for the relativistic electrons in the ground of quantum hydrodynamics model and the 
results are compared to the case of semiclassical Thomas-Fermi dusty plasma. Using 
the standard pseudo-potential approach it is remarked that the Fermi-Dirac plasma, 
in contrast to the Thomas- Fermi counterpart, accommodates a wide variety of non- 
linear excitations such as positive-/negative-potential ion solitary and periodic waves, 
double-layers and double-wells. It is also remarked that, the relativistic degeneracy 
parameter which relates to the mass-density of plasma has significant effects on the 
allowed matching-speed range in dusty Fermi-Dirac plasmas. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. BACKGROUND 



The study of plasma under extreme conditions, such as ultrahigh pressure or temperature 
is growing rapidly and becoming one of the intense fields of plasma research. There are 
already numerous laboratory applications for plasmas under high pressure such as that 
found in semiconductors and nanomaterials l|. On the other hand, plasmas under high 
temperature are the main ingredients of inertial confinement fusion and Tokamaks 2j-y|. 
It is known that many astrophysical places are under such extreme conditions The 
study of ultradense plasmas, however, requires the relativistic and quantum mechanical 
considerations. For instance, the quantum effects become more important when the Pauli 
exclusion principal acts effectively on fermions and the interparticle distances get much lower 
than the characteristic de Broglie thermal wavelength, /i/(27rm/csT) 1 / 2 [11] which leads to 
the degeneracy condition. Under a very high external pressure, for instance, the degenerated 



electrons in matter become relativistic and obey the more general Fermi-Dirac statistics [12 ]. 

It has shown by S. Chandrasekhar [13] that in a white dwarfs the relativistic pressure, 
under the influence of gravity, may eventually lead to the collapse of star, due to softening 
of degeneracy pressure. Using the Fermi-Dirac statistics for electrons, a rigorous mathe- 
matical criteria has been developed 14| under which a white dwarf can be considered as 
completely degenerate and has been shown that a typical white dwarf can be examined suc- 
cessfully under the ideal zero-temperature Fermi-gas model such as that observed in ordinary 
metals. New study based on semiclassical Thomas-Fermi model reveals that the nonlinear 
ion-acoustic wave dynamics such a their propagations and collisions in superdense magne- 
tized plasmas are much different when the electrons are ultra-relativistically degenerate from 
those with non-relativistic electrons 
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16j. Therefore, it would be interesting to explore 



the effects of relativistic degeneracy degree (relativistic degeneracy parameter) on nonlinear 
wave dynamics in a white dwarf under high gravitational pressure. 

There has been already many investigations on the properties of nonlinear propagations 
in a degenerate plasma which examine either non-relativistic U2I] or ultra-relativistic 
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26| degeneracy conditions. However, the later studies in ultra-relativistic degeneracy 



cases are based on a much simplified model of Thomas-Fermi approximation which has its 
known limitations. In current work, using the quantum hydrodynamics formulation [27| and 
assuming a more appropriate Fermi-Dirac model of degeneracy pressure which accounts for 
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the relativistic effects of electrons, we examine the general case which covers the whole range 
of the relativistic degeneracy, special limits of which should produce the previous results. 
The presentation of the article is as follows. The basic normalized plasma equations are 
introduced in section [Til Nonlinear arbitrary-amplitude solutions are derived in section IIHI 
Numerical analysis and discussions is given in section IIVI and final remarks are presented in 
section [V] 



II. BASIC FLUID MODEL 



In order to evaluate the nonlinear dynamics of ion waves we use the complete set of 
quantum hydrodynamics (QHD) equations. Here, we consider a dense plasma consisting 
of relativistically degenerate cold electrons, dynamic ultra-cold ions and neutralizing static 
background dust ingredients. It is assumed that such a plasma is collision-less due to Pauli- 
blocking. It is further required that the electrons obeythe zero-temperature Fermi-Dirac 



statistics. Therefore, the closed set of QHD equation 28j, in dimensional form, is written 
as following 



drij _|_ dnjVj _ n. 

dt " r dx U ' 

dn e I dn e v e n 

dt " r dx U ' 



dvj _i_ y.dvi §_dj>_ _ ^k B Tj I _n^_ \ drn / j\ 

dt 1 dx rrii dx m; n^o \ n i0 J dx ' v ' 

UhL (9Ve _L_ y §Ve\ _ _e_ d± 1 8P e , h 2 d ( 1 9 2 y/T 

rrii \ dt e dx I rrii dx min e dx 2m e rrii dx \ yfn2 dx' 2 

= 4ire(n e - m + z d n d ), 

where, and v j denote the j-th species number-density, mass, atomic-number and 

velocities with j = {e,i,d} for electron, ion, dust, respectively. Note that the quantity v e 
stands for the relativistic electron speed and n^o denotes the ion equilibrium density. Also 
7, Tj and h, respectively, present the adiabatic constant, ion-temperature and the scaled 
Plank-constant. Furthermore, the following scalings can be used in order to obtain the 
normalized set of equations 



x ->■ —x, t -> — , n e ,i -> n eji n , v e4 -> v £ti c s , -> 4>'-^—, (2) 



c s t -m e c 



where, the bar symbols denote the dimensionless quantities which are omitted in forthcom- 
ing algebra for simplicity. The parameters u P i = \jATie 2 nio/mi and c s = \J m e c 2 / rrii are 
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characteristic plasma frequency and sound-speed, respectively, and c is the speed of light in 
vacuum. Also, no denotes the normalizing number- density, no = &n ^f — 5.9 x 10 29 cm~ 3 
a value in white dwarf electron number-density range. The electron relativistic degeneracy 
pressure in zero-temperature configuration, which governs most dense astrophysical objects 
such as white dwarfs, can be expressed in the following general form [13 ] 

4 5 

TTTfi C ( r ~1 1 

P e = _^|i?(2i? 2 -3)v / rT^ + 31n[ J R + v / rT^J), (3) 

in which R = [n e /no) 1 ^ (no = 8n ^3 C — 5.9 x 10 29 cm~ 3 ) is the relativistic degeneracy 
parameter measuring the relative Fermi- momentum of electrons, i.e. R = pp e /m e c. In the 
limits of very small and very large values of the relativity parameter, R,we obtain the well 
known non- and ultra- relativistic degeneracy pressures (e.g. see Ref. 13J) 

P e = { 20 H m % e v " }. (4) 




Using the scalings introduced above, the normalized hydrodynamics set of basic equations, 
including the quantum force effect, is written as 

at " r dx u ' 

dn e I dn e v e n 

dt " r dx U ' 

dvj I v .§Vi_ _ _ ( lk B Tj \ ( ni_\' 1 dni /g\ 

dt 1 dx dx \ riiQin^c 2 J \ mo / dx ' * ' 

1 d 2 y/n, 



a 2 <a 

^£ = n e - 7i» + ^n rf , 

where, i?o = ( n io/ n oY^ 3 is a measure of the relativistic degeneracy effects (normalized relativ- 
ity parameter), a = n e o/n i0 and (m e /mj)i/ 2 /2 = ft 2 o; 2 j/(m e c 2 ) 2 , is the quantum diffraction 
coefficient which is related to the relativity parameter via (m e /mj)iJ 2 /2 ~ 4.6 x 10~ 10 Rq. It 
is clearly confirmed that the quantum correction in this model is negligible even for higher 
values of relativity parameter. For instance, for the value of Rq = 1000 we obtain a quantum 
correction of order ~ 0.46. It should also be noted that the relativity parameter, Rq, (using 
p ~ 2m p (nio + ZdUd) = 2m p nio(2 — a) with m p the proton mass) is related to the mass density 
through p ~ 3.94 x 10 6 (2 — a)R^(gr / cm?). Taking, for instance, the density of typical white 
dwarfs to be in the range 10 5 < p(gr/cm 3 ) < 10 8 , results in typical values of 0.37 < Rq < 8 
for the relativity parameter. However, the current study does not limit the applicability of 
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results to specific density range. Now, neglecting the small contributing terms containing 
the ratio m e /mi and assuming that ksTi <C m e c 2 , we are led to the following simplified set 
of normalized equations 

drij _|_ dnjVj _ rv 

at " r dx u ' 

dn e _|_ i9n e n e _ n 
8t " r 9x u ' 

at ^ ^ dx dx ' v ; 

g - IVi + *§« 2/ V/ 3 = o, 

a 2 4> 



^ _ n e - m + (1 - a) 



Note, however, that the overall neutrality condition, (zdndo)/n i0 = (1 — a), is met at <fi = 0. 
Solving the basic equations (Eqs. (jSJ)) for the electron and ion number-densities in terms of 
electrostatic-potential, using the quasi-neutrality condition, and employing the appropriate 
boundary requirement (limrt e j = {a, 1}), results in the following energy relations 

ILi r, w 

« 1 ~m* (7) 



n e = a 1 R 



i 3/2 



R 2 a 4 ^ 3 + (2^/1 + i?g a 4/3 + 

The energy relation for electrons resembles to that of known Thomas-Fermi 
semiclassical model (n e = a(l+(p) 3 ^ 2 ) used in some literature jl5l.ll61.l23l. 29]. How- 



ever, this new expression contains more information about the electron dynam- 
ics as it contains an extra parameter which describes the relativistic dynamical 
properties of electrons (electron relativistic Fermi-momentum , i.e. R = pF e /m e c) 
in an ultra-dense plasma. It should be noted that the Thomas-Fermi approxima- 
tion is valid only for slow variation of electrostatic potential and its application 
is limited for the ultradense plasmas where electron exchange potential can be 
dominant and one should appeal for the Thomas-Fermi-Dirac approximation 
{24}]. An advantage of this model with respect to Thomas-Fermi is that the cur- 
rent model accounts for the whole range of relativity parameter the two limits 
of which leads to the well known non- and ultra relativistic degeneracy pressure 
limits given in Eq. 

In the proceeding section we aim at finding the appropriate pseudo-potential 
which describes the propagation of nonlinear ion-acoustic waves in the plasma 
described by Eq. ([6]) and compare the results with that obtained in the Thomas- 
Fermi model. 
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III. ARBITRARY-AMPLITUDE ANALYSIS 



In this section we aim at finding the appropriate Sagdeev potential in plasma defined 
by Eq. (jSj) which satisfies the nonlinear ion-solitary (ISWs) as well as periodic (IPWs) 
solutions. To do this using the standard procedure, we transform into a moving coordinate 
relative to the wave-train. Thus, using the new coordinate £ = x — Mt in order 
to find the stationary solution with M being the matching-speed of stationary 
frame not to be confused by the Mach number (e.g. see Ref. j^sj]), a measure 
of the wave-speed relative to the normalizing parameter, c s , and substituting Eqs. 
(J7J) into Poisson's equation, multiplying by and integrating once with suitable boundary 
conditions mentioned in the previous section, we derive the well-known energy integral 

2 



2'lde 



0. 



(8) 



where, the desired Sagdeev potentials U(<f)) for the Fermi-Dirac model, read as 
U FD {4>) = -8aR% \(1 - a)<j> + My/M 2 - 2^1 + R a 2 / 3 (2R 2 a 4 / 3 - 3) y/l + Rfc 4 ' 3 



3 (2(f) 2 - 1) y/l + R 2 a 4 / 3 + <P (2<p 2 + 1) + 2R 2 a A / 3 (30 + y/l + R 2 a 4 / 3 



x ^R 2 a^ 3 + <p (V + 2^1 + R 2 a A /^j + Ssinh" 1 (i? « 2/3 ) + 3 In 2 + 8aM 2 R t 
-3 In 



(9) 



2 (V + y/l + R 2 a^ + 2jR 2 a^ + (<j> + 2y/l + R 2 a A l 3 



Similarly, for the Thomas-Fermi plasma of this kind, we obtain 



r 77- i<>) = M 2 + ^ - <j) + a<j) - |a(l + <f)) 5/2 - M 2 v/l 



(10) 



5 ' T ' "~ V " M 2 ' 

It is remarked that the values of U(4>) and <9^?7(0) both vanish at <p = for both 
pseudo-potentials, as expected. From Eqs. ([7]) we may also derive the maximum 
values of <p in terms of other plasma parameters for Fermi-Dirac case, as 



%i 



M 2 
2 ' 



(11) 



m2 = ±l-y / l + J R 2 a 4 / 3 , (12) 
and one maximum and minimum values of <p for the Thomas-Fermi case, as 



M 2 
~2 



-1, 



(13) 



The possibility of solitary waves then requires the following conditions to met, simultane- 
ously 



dU{4>) 



0=0 



0. 



d 2 U{ct)) 



4>=o 



< 0. 



(14) 



0=0 



It is further required that for at least one either maximum or minimum nonzero 0-value, 
we have U(<f> m ) = 0, so that for every value of in the range m > > or m < < 0, 
U(4>) is negative. In such a condition we can obtain a potential minimum which describes 
the possibility of ISWs propagation. 

On the other hand, the existence of a nonlinear periodic wave requires that 



dU(cf>) 



0=0 



0. 



d 2 U((j)) 



0=0 



> o, 



(15) 



0=0 



It is also required that for at least one either maximum or minimum nonzero 0-value, we 
have U(4> m ) = 0, so that for every value of in the range m > > or m < < 0, 
£7(0) is positive. In such a condition we can obtain a potential maximum which describes 
the possibility of IPWs propagation. 

Since, the Sagdeev potential has triple roots, there is another possibility of obtaining 
double-layer (DLs) or double- wells (DWs). The condition for a solitary (periodic) DLs 
is W)l0=0 m = ^(0)/#|^ m = and d 2 U {4>) / d4>\ =(f>m < {d 2 U {4>) / d4>\ = ^ > 0). 
Existence of a DW, on the other hand, requires that £7(0 m i) = U((p m 2) = 0, where, m i and 
m 2 are oppositely signed. A DW is condition in which rarefactive (0 m < 0) and compressive 
(0 m > 0) nonlinear excitations coexist. Since, evaluation of the complete set of conditions 
on the pseudo-potential analytically is not possible in current case, in the proceeding section 
a numerical scheme will be employed. 



IV. DISCUSSION 

Detailed numerical investigation of the Sagdeev potential Eq. ([9]) indicates that a wide 
variety of nonlinear excitations can be present in a dusty Fermi-Dirac plasma define in Eq. 

Figure 1 shows positive and negative potential arbitrary amplitude ISWs and IPWs and 
the variation of the potential profiles with respect to relative electron concentration, a, for 
fixed values of matching-speed, M, and relativity parameter, R (or plasma mass-density). It 
is remarked that the increase/decrease in fractional electron concentration (decrease/increase 
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in relative dust concentration) for rarefactive ISWs leads to increase/decrease in soliton 
amplitude/width, while for IPWs this increase has only minor effects (increase/decreases) 
on the width/amplitude of the soliton. For IPWs, on the other hand, the increase of relative 
electron number-density leads to the decrease/increase of wave amplitude/ width for both 
positive and negative potential (<f> m < 0) waves. 

Figure 2 shows the formation of W-shaped positive and negative potential periodic-DLs 
(Figs. 2(a) and 2(b)) and solitary-DWs (Figs. 2(c) and 2(d)) in a Fermi-Dirac dusty plasma. 
From Fig. 2(c) it is clearly observed that the compressive soliton shapes are affected via 
the variation of the matching-speed more than that of rarefactive ones. On the other hand, 
Fig. 2(d) reveals that the variation of relative electron concentration effects the rarefactive 
solitary excitation profiles more than that for compressive ones. It is also noted that the 
increase of a on rarefactive and compressive solitary profiles, discussed in previous paragraph 
are quite similar to those of solitary-DWs. 

In Fig. 3 we showed the variation allowed matching-speed range for ISWs 
with respect to the normalized plasma mass-density, p, for the case of Fermi- 
Dirac plasma (using the relation p = p/p ~ (2 — a)R^ with p = 3.94 x 10 6 (gr/cm 2 )) 
which is obviously not possible for the case of the Thomas-Fermi plasma. It is 
remarked that, the speed range of possible propagation of ISWs increases as the 
normalized plasma mass-density increases for all values of the relative electron 
concentration. Note that the numerical values of c s is much higher (compared to 
the electron Fermi-Speed) in here compared to that in ordinary non-degenerate 
plasmas. It is also revealed from Figs. 3(a)-3(d) that the increases of the relative 
electron concentration in the Fermi-Dirac plasma leads to the increases in the 
allowed range of the matching speed. 

The regions of possibility of these nonlinear excitation are shown in colors in Fig. 4 with 
each color representing a specific type of nonlinear excitation in the Fermi-Dirac plasma. 
There is also a plot (Fig 4(d)) with numerated regions for non-color option. The description 
of different region for nonlinear excitations is as follows 

Region- 1 (in blue color) indicates the possibility of positive potential IPWs an example of 
which is shown in Fig. 1(c). 

Region-2 (in black color) indicates the possibility of solitary-DWs an example of which is 
shown in Fig. 2(c) and 2(d). 
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Region-3 (in red color) indicates the possibility of negative potential IPWs an example of 
which is shown in Fig. 1(d). 

Region-4 (in yellow color) indicates the possibility of compressive ISWs an example of 
which is shown in Fig. 1(b). 

Region-3 (in green color) indicates the possibility of rarefactive ISWs an example of which 
is shown in Fig. 1(a). 

The periodic-Dls, examples of which is shown in Figs. 2(a) and 2(b), exist at the left 
borders of regions 1 (yellow), 2 (black) and 3 (red), where, at the left borders of regions 
1 (yellow) and 2 (black) there exist positive potential-DLs and the negative potential-DLs 
occur at the left border of region 3 (red). It is further remarked from Fig. 3 that, there is 
a region for approximately 3.75 > a > 0.05 for which no IPWs are present and that this 
range is affected only slightly with the variation of the relativistic degeneracy parameter. 
However, as the parameter R which represents the mass-density of the plasma increases 
the regions expand nearly in uniform manner and move to higher Mach-number values. 

Finally, in Fig. 5, we show the regions of existence for ISWs (region in yellow 
labeled 4) and IPWs (region in red labeled 3) for the corresponding dusty plasma 
in Thomas-Fermi model. It is clearly remarked that, only two regions, namely 
compressive ISWs and negative potential IPWs regions, are present in the case 
Thomas-Fermi dusty plasma. 

V. CONCLUDING REMARKS 

The Sagdeev pseudo-potential approach was used to investigate the propagation of non- 
linear ion excitations in a dusty Fermi-Dirac electron-ion plasma with considering the effect 
of electron relativistic degeneracy. The matching criteria of existence of such solitary excita- 
tions were numerically investigated in terms of the relativity parameter, R , the relative dust 
concentration, (1 — a), and the allowed matching-speed, M and there was found five dis- 
tinct regions which accommodate a variety of nonlinear excitations. Nonlinear double-wells 
and double-layers are also found to exist in such plasma. Furthermore, it was confirmed 
that the relativistic degeneracy parameter has significant effect on the range of allowed 
matching-speed in dusty Fermi-Dirac plasmas. The current findings can be applicable to 
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highly degenerate dense astrophysical compact objects such as white dwarfs. 
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FIGURE CAPTIONS 



Figure- 1 

(Color online) Seudo-potential profiles in Fermi-Dirac dusty plasma and their variations 
with respect to the fractional electron concentration, a, for fixed other plasma parameters. 
The dash size is correspondingly varied with the change in the value of a parameter. 

Figure-2 

(Color online) The formation of positive-/negative-potential double-layers (DLs) and 
double-wells (DWs) in Fermi-Dirac dusty plasma. The dash size is correspondingly var- 
ied with the change in the value of a or M parameters. 

Figure-3 

(Color online) The variations of allowed ISW matching-speed with respect to normalized 
plasma mass-density, p, and the relative electron concentration, a in the Fermi-Dirac dusty 
plasma. 

Figure-4 

(Color online) Regions of existence of different nonlinear excitations shown in colors 
(blue-black- red- yellow-green) enumerated, respectively, (1-2-3-4) in Fig. 3(d), in Fermi- 
Dirac dusty plasma. Each region belong to a distinct kind of nonlinear wave (refer to text). 
Variations in the regions with respect to the relativistic degeneracy parameter, R , which 
also reflects the plasma mass-density through the relation, p ~ 3.97 x 10 6 (2 — a)Rl(gr/cc), 
is shown through plots (a)-(d). 

Figure-5 

(Color online) (a) Regions of existence for ISWs (region in yellow labeled 4) and IPWs 
(region in red labeled 3) and (b) the variation of pseudo-potential for ISWs with respect 
to the change in relative electron concentration for the fixed matching-speed, in a Thomas- 
Fermi dusty plasma. 
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